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We review the models developed and techniques used 
in recent years to study the kinetics of phase ordering in 
a class of complex fluids, namely, ternary microemulsions 
and micellar solutions. 



I. INTRODUCTION 

The "head" part of surfactant molecules consist of a 
polar or ionic group. The "tail" of many surfactants con- 
sist of a single hydrocarbon chain whereas that of some 
other surfactants, e.g., phospholipids, are made of two 
hydrocarbon chains both of which are connected to the 
same heada. In contrast, gemini surfactantaSlj, consist 
of two single-chain surfactants whose heads are connected 
by a "spacer" chain and, hence, these "double-headed" 
surfactants ate sometimes also referred to as "dimeric 
surfactants"BLl (see fig.l). 




igle-chain Double-chain 



hydrophilic head 



hydrophobic tail 



taining water. These not only find wide ranging appli- 
cations in detergent and pharmaceutical industries, food 
technology, petroleum recovery, etc. but are also one of 
the most important constituents of cells in living sys- 
tems. Therefore, physics, chemistry, biology and tech- 
nology meet at at the frontier area of interdisciplinary 
research on association colloids formed by surfactantsllil. 

A ternary microemulsion is a "colloidal" complex fluid 
consisting of three components, namely, oil, water and 
surfactants; the phase diagrama-jjif such ternary systems 
have been studied extensivelyHEl. When the concen- 
tration of the amphiphilic molecules is sufficiently high, 
the low-temperature equilibrium phase of the system is 
lamellar where the amphiphiles arrange themselves in 
(approximately) parallel stacks. On the other hand, 
when the concentration of the amphiphiles is not high, 
the system exhibits either a droplet phase or a bicontinu- 
ous microemulsion phase in equilibrium at sufficiently low 
temperatures depending on the relative concentrations of 
4 water. If the concentrations of oil and water are 
able then the bicontinuous microemulsion phase is 
ed. But, if the concentration of oil (water) is much 
than that of water (oil) then droplets of oil (water) 
bund to be dispersed in water (oil); these droplets 
ften referred to as micelles and, therefore, the sys- 
under such conditions are called micellar solution, 
se systems find wide ranging industrial applications, 
e aim of this chapter is to present a systematic and 
Gemirii'^^^ date review of the models developed and techniques 
used so far to study the kinetics of phase ordering in-Jiii-. 
croemulsions and micellar solutions (see referencesll3"t3 
for earlier reviews). 

spacer 




FIG. 1. Different types of amphiphiles. 



II. KINETICS OF ORDERING IN BINARY 
MIXTURES IN THE ABSENCE OF 
SURFACTANT IMPURITIES 



When put into an aqueous medium, the "heads" of the 
surfactants like to get immersed in water and, hence, are 
called "hydrophilic" while the tails tend to minimize co 
tact with water and, hence, are called "hydrophobic" 
The spacer in gemini surfactants is usually hydropho- 
bic but gemini surfactants with hydrophilic spacers have 
also been synthesizedB. Surfactant molecules are called 
"amphiphilic" because their heads are "water-loving" 
and hydrocarbon chains are "water-hating" . Because of 
their amphiphilicity the surfactant molecules form ^'self- 
assemblies" (i.e., supra- molecular aggregates), such as 
monolayer and bilayer membranes, micelles, inverted- 
micelles, etcu, in a multi-component fluid mixture con- 



We introduce the key concepts involved in the studies 
of phase ordering by illustrating them with the help of bi- 
nary mixtures. The techniques used so far for the studies 
of phase ordering in ternary microemulsions and micel- 
lar solutions are extensions of those used widely for the 
corresponding studies in the simpler case of binary mix- 
ture of immiscible components, e.g., binary alloys, binary 
fluids, etc. Therefore in this section we classify and. t 
marize the main techniques used for such studiealj Ei 

Suppose a binary alloy has been quenched from a high 
temperature Th to a temperature Ti well below the coex- 
istence curve. The morphology of the coarsening domain 
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"pattern" depends on the relative concentrations of the 
two components. If the concentrations of the two com- 
ponents are comparable then the pattern has a random 
"interconnected" structure. On the other hand, if the 
concentration of one of the components is much smaller 
than that of the other the pattern consists of "droplets" 
of various sizes. 

There are essentially two different (albeit complemen- 
tary) theoretical approaches to the study of the kinetics 
of the phase ordering processes, viz., microscopic models 
and phenomenological models. 

A. Microscopic Models 

Two types of microscopic models have been used; the 
molecular models defined on a continuum and lattice 
models. 



1. Molecular Models on a Continuum: 

For the Molecular Dynamics simulation of an immisci- 
ble binary mixture of two species, labelled A and B, one 
first postulates an approximate form of inter-molecular 
interactions, e.g., truncated Lennard- Jones potentials. 
The potentials are chosen in such a manner that the in- 
teractions A — A and B — pS are attractive while the inter- 
action A — B is repulsivcEHl. Normally one uses the Verlet 
algorithm for integrating Newton's equation of motioned. 
In principle, one can use either a constant-energy ensem- 
ble or a constant temperature ensemble. 

2. Microscopic Lattice Models 

The binary alloys can be modelled as an Ising spin 
system where Si = 1 correspond to an A atom and Si — 
— 1 correspond to an B atom. In the symmetric case, 
i.e., when Eaa = Ebb, the Hamiltonian for the system 
is given by 

H = -JY.S,S, (1) 

where J > and the summation is to be carried out over 
all the distinct nearest-neighbour spin pairs. However, 
no such model is complete without specification of the 
prescription for the dynamical evolution of the systemEJ. 
In the case of binary alloys, unlike the magnetic coun- 
terpart, the order parameter is conserved, i.e., the con- 
centrations of A atoms and B atoms remains unaltered 
during the time evolution in a closed system. The ki- 
netics of ordering in such a system can be studied at a 
microscopic level by using the so-called Kawasaki spin- 
exchange dynamics: two anti-parallel nearest-neighbour 
spins can exchange their position with a probability 
1/[1 + exp(/3A£;)], with /3 = l/lfcsT) where fcs is the 



Boltzmann constant and A_B is the energy change that 
would be caused by the interchange of the two spins. 

B. Phenomenological Models with Langevin 
Dynamics 

In the phenomenological approach the system is de- 
scribed by an order parameter field '(/'(r), which is the 
local difference in the concentrations of the A and B 
atoms. In contrast to the discrete allowed values of the 
spins in the lattice model, the order parameter can take 
all real values in the interval — 1 < '0 ^ 1- The coarse- 
grained free-energy functional for the d-dimensional sys- 
tem is given by 

F[i/;(f)] = j d'*r[roV'^(r) + uip^{r) + c\Vi}{f)\^] (2) 

where tq, u and c are phenomenological coefficients. The 
symmetry requirements rule out the possibility of tp and 
tp^ terms in this functional. The dynamics of the sys- 
tem is assumed to be governed by the so-called Langevin 
equation 

dij{r, t) /dt ^ rV^ [5F/5^{r, t)] -f 7^{r, t) (3) 

where F is the phenomenological kinetic coefficient, 5/ Sip 
denotes the functional derivative with respect to ■0 and 
the Laplacian takes care of the fact that the order pa- 
rameter is conserved, rjij^ t) is the noise, which is usually 
assumed to be of 'Gaussian white' nature, i.e., 

< ?7(f, t)r^{r', t') >= 2kBTV^S{r- r')5{t - t') (4) 

where ksT guarantees the approach to the true Gibb- 
sian equilibrium. In the context of the binary alloys the 
equation (3) is called the Cahn-Hilliard-Cook equation 
for historical reasons. The model described by the.equa- 
tions (2) - (4) is usually referred to as the model SEj. The 
computation becomes much more efficient (i.e., the mor- 
phological characteristic of the asymptotic regime can be 
obtained with a very short computer time) byj-solving 
equation (3) above by the cell dynamics methodcJ. 

C. Characteristic Quantities of Interest 

Two quantities are most important in describing the 
kinetics of this growth process, viz., the time dependence 
of the typical linear size R{t) of the ordered regions and 
the dynamic scaling of the structure factor S{q,t) char- 
acterizing the statistical self-similarity of the coarsening 
pattern at different times. The structure factor is the 
Fourier transform of the two-point correlation function 
in real space, i.e., 

S{q, t)^J2 ^(^^' t)^Mi'l-r1 (5) 
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where G{f,t) =< S{6,t)S{f,t) > - < S is the cor- 
relation function in real space. However, since numerical 
computations are carried out either on discrete lattices 
(in the case of microscopic models) or on discrete grids (in 
the case of phcnomcnological models) the corresponding 
structure factor is given by 

S{q,t) =< \{l/N)Y,G{n,t)expiiq.n)\^ > (6) 

with q = {2tt/ L)(mx + ny+pz) and m,n,p = 1, 2, 3, L. 

Following a rapid quench from a very high temperature 
to a temperature below the coexistence curve, coarsening 
of ordered domains takes place and, consequently, the 
first zero crossing of G(i?, t) (i.e., the smallest R for which 
G{R, t) = 0) occurs at larger and larger values of R at 
successively longer values of time t, where G{R, t) is the 
circularly averaged two-point correlation function. As a 
result, the location of the first zero crossing of G{R, t) 
may be taken as a measure of R{t). 

During the coarsening process the position of the peak 
in S{q, t) keeps moving towards smaller values of q where 
S{q, t) is circularly averaged. The dynamical scaling form 
of the structure factor is given by 

S{q,t)^R''F{qR{t)) (7) 

where F{x) is a function of x = qR{t). The length R{t) 
can be extracted from S{q,t) in several different ways: 
(i) R~^{t) — A:„j/27r, where km is the location of the 
maximum of the structure factor S{q,t), (ii) 2it/R may 
be identified with the first moment of the structure factor 
or the square root of the second moment of the structure 
factor, etc. These measures of R{t) can also be used in 
the case of ternary microemulsions and micellar solutions. 
Moreover, one can compute the mean-cluster size of a 
particular component, say A, using the definition 

XA^[J2 sMmT.'^('^ (8) 

where n(s) is the number of clusters of type A with s sites 
and Smax is the corresponding size of the largest cluster. 

For binary alloys, for which the order parameter is con- 
served, i.e., for model B, the length R{t) has been found 
to follow the growth law 

R{t) - r (9) 

where n = 1/3. In contrast, for systems with non- 
conserved scalar order parameter (model A) n ~ 1/2. 
Moreover, in case of binary fluids hydrodynamic effects 
lead to the growth law (9) with n = 1. Furthermore, 

R{t) ^ {logtr, (10) 

with dimensionality-dependent exponent x, in phase sep- 
arating binary systems-in the presence of quenched dis- 
order (e.g., impurities)EZI. 



III. KINETICS OF ORDERING IN TERNARY 
MICROEMULSIONS AND MICELLAR 
SOLUTIONS 

In our discussion of the kinetics of phase ordering in 
ternary microemulsions and micellar solutions we shall 
begin with the molecular models, which are truly micro- 
scopic description of the system, and then consider mod- 
els at a coarser level and, finally, consider phcnomcnolog- 
ical models. 



A. Molecular Dynamics Simulations of Molecular 
Models on a Continuum 

Laradji et a/.El have carried out a MD simulation of 
phase separation in a binary mixture in the presence 
of surfactants. The interaction potential between two 
molecules i and j (oil-oil, water-water or oil-water) is as- 
sumed to be 

f/(rl,) = 46[((7/r,,y2 _ (2^^^^^ „ l)(<T/r,,)'] (H) 

where 6aiaj is the Kronecker delta function and ai de- 
notes the type of the molecule i, i.e., ai ~ 1 for a water 
molecule and 2 for an oil molecule. Thus, the potential 
above ensures that molecules of different species always 
interact repulsively. For example, if ai ^ aj, 

Uin,) = Mia/r,,y^ + ia/n,f] (12) 

On the other hand, if ai — aj, 

U{n,) = 4e\ia/n,r-ia/n,f] (13) 

In this model the "water-loving" head of every surfac- 
tant molecule is treated as identical to a water molecule 
and the "oil-loving" tail is treated as identical to an oil 
molecule. In other words, every surfactant molecule is 
assumed to be, effectively, a diatomic molecule one part 
of which is water-like and the other part is oil-like; the 
two parts of the molecule are assumed to be connected 
to each other by a harmonic spring so that the potential 
is given by 

f/,,(f) = (X,/2)(r„--4)'. (14) 

To our knowledge, this type of models W|exe first proposed 
by Smit et al.tB. However, Laradji et al!E3 was the first to 
study the kinetics of phase ordering by MD simulation of 
such a model. A non-zero concentration of amphiphiles 
was found to slow down the coarsening process leading 
to a significant deviation from the power-law (9). 

B. Microscopic Lattice Models 
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1. Widom Model: a Spin-1/2 Ising Model on Simple-Cubic 
Lattice with Farther-Neighbour Interactions 

Widom modeS is the simplest lattice model of ternary 
mixtures of oil, water and surfactants. This is a lat- 
tice model in the same spirit as, for example, the lat- 
tice model of binary alloys. However, in contrast to the 
lattice model of binary alloys described in the preced- 
ing section, the molecules of oil, water and surfactants 
in this model are located on the bonds rather than on 
the lattice sites. In this model the molecules of oil, water 
and surfactants are represented by the nearest-neighbour 
bonds of a spin-1/2 Ising model where every spin in- 
teracts with the nearest-neighbour as well as a specific 
subset of farther-neighbour spins. The bonds between 
the up-up nearest-neighbour spin pairs represent water 
molecules, those between down-down nearest-neighbour 
spin pairs represent oil molecules and those between an- 
tiparallel nearest-neighbour spin pairs denote the surfac- 
tant molecules. The Hamiltonian for this model is given 
by 

H = -JY^ SiSj - 2M S,Sk ~ M S^Si (15) 

<ij> <ik> <il> 

where, for three dimensional systems, the summations in 
the first, second and third terms on the right hand side 
are to be carried out over the nearest-neighbour, second- 
neighbour and fourth-neighbour spin pairs, respectively, 
on a simple cubic lattice. The interaction J is positive 
(ferromagnetic in the terminology of spin models of mag- 
netic materials) whereas M is negative (antiferromag- 
netic). The farther-neighbour intejLactions arise from a 
prescription, suggested by Widon£3 , to take into ac- 
count the bending rigidity of the amphiphilic monolayer 
membrane at the oil-water interface. This model reduces 
to the standard spin-1/2 Ising model with only nearest- 
neighbour ferromagnetic interactions on a simple cubic 
lattice if A/ = 0. 

For studying r-the dynamics of the Widom model 
Morawietz et a/.til introduced a Kawasaki-type spin ex- 
change dynamics which conserves the molecules of oil, 
water and surfactants. Note that in this algorithm the 
spins on the lattice sites, rather than the molecules 
on bonds, are exchanged; however, the algorithm is 
such that it leads to simultaneous exchange of three 
molecules while satisfying the condition that the number 
of molecules of each species remains conserved. The al- 
gorithm is as follows: a up-spin can exchange its position 
with one of its down-spin neighbours, with probability 
1/(1 -f exp(PE)), provided the numbers of up and down 
spins among the neighbours of the first spin (excluding 
the second spin, which is also a neighbour of the first) are 
the same as the neighbours of the second spin (excluding 
the first spin, which is also a neighbour of thCpSecond). 

Using this new algorithm, Morawietz et alEj investi- 
gated the kinetics of dis-ordering, rather than the tradi- 
tional study of ordering, by reverse quenching the system 



from the ordered phase to the disordered phase. More 
specifically, an initial configuration was created where a 
layer of surfactants in the central part of the system is 
sandwiched between a column of oil above and a column 
of water below. The parameters J/{kBT) and M/{kBT) 
were so chosen that the corresponding equilibrium phase 
is known to be a disordered fluid. To our knowledge, the 
kinetics of ordering in the Widom model has never been 
studied. 



2. K aw akatsu- Kawasaki Model: a Decorated Spin-1/2 Ising 
Model 

The Alexander modelsl is a decnrated spin-1/2 model, 
which was extended by Chen et a/H for a more realistic 
descripticftL. of ternary microemulsions. Kawakatsu and 
KawasakiEfl have studied the kinetics of phase ordering 
using a decorated spin-1/2 model which may be regarded 
as a simplified special case of the model considered by 
Chen et a/Ej. In this model every lattice site is occupied 
by a Ising spin; Si = 1 corresponds to a water molecule 
and Si — —1 corresponds to a oil molecule at the i-th 
site. Just as in the molecular models mentioned in the 
preceding section, the "water-loving" head part and the 
"oil-loving" tail part of each surfactant is assumed to be 
identical to a water molecule and an oil molecule, respec- 
tively. A fraction of the nearest-neighbour pairs of dis- 
similar molecules are connected by a rigid bond; each of 
the dumb bell-like structures where a rigid bond connects 
two dissimilar molecules at its two ends is identified as a 
surfactant. Thus, in this model, the molecules of oil and 
water occupy the lattice sites whereas the surfactants are 
located on the bonds. The Hamiltonian for this model is 
given by 

H= (J/2) ^(1-5.5,) (16) 
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pairs considered for exchange are next-nearest neighbours 
with respect to each other (see fig. 2 for the details). 
C^^¥t# most important observation of this Monte Carlo 
study is that the surfactants slow down the growth pro- 
cess; this is consistent with the results of the MD simu- 
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latio^ mentioned in the preceding section. 



FIG. 2. Schematic representations of (a) the lattice model of 
a binary alloy and (b) the Kawakatsu- Kawasaki lattice model of 
a ternary system where one of the components is a surfactant. 
The elementary 'exchange' processes allowed in the dynamics of 
the two models are also shown. Exchange of anti-parallel nearest 
neighbours only is allowed in the case of the binary alloy whereas 
exchange of antiparallel next-nearest neighbours only is allowed 
in the case of ternary system 



3. Larson Model: Self-Avoiding Chain Model of Amphiphiles 

^ If© Larson modelH^l was originally developed for 
ternary microemulsions. The fluid under investigation 
is modelled as a simple cubic lattice of size x Ly x Lz- 
Each of the molecules of water (and oil) can occupy a 
single lattice site. A surfactant occupies several lattice 

sites each successive pairs of which are connected by a 

rigid nearest-neighbour bond. A^single-chain surfactant 
cm be described by the symbolt3 %nJ\fp7iq where T de- 
nptes tail, Ti denotes head and TV denotes the 'liaison' 
or neutral part of the surfactants, to, p and q are inte- 
gers denoting the lengths of the tail, neutral region and 
head, respectively, in the units of lattice sites. Thus, each 
^pgle-chain surfactant is a self-avoiding chain of length 
I = {m + p + q). We shall refer to each site on the sur- 

f ptants as a monomer. The "water-loving" head group 

isassumed to be "water-like" and, similarly, the "oil- 
loving" tail group is assumed to be "oil-like" . Figure 3 
shows Larson-type models for different types of surfac- 
d mts. 



Note that the energy of the water-water, oil-oil and 
water-oil interactions are given by fi'mtu = = Eoo and 
Ewo — J- Since J > 0, water and oil tend to phase sepa- 
rate in the absence of surfactants. Since the two ends of 
the surfactants are not allowed to split into two separate 
water and oil molecules during the time evolution of the 
system the interaction energy between the water-like end 
and oil-like end of such a surfactant does not affect the 
dynamics. However, the dynamics must distinguish be- 
tween individual water (and oil) molecules and the water- 
like (and oil-like) ends of surfactants so as not to split a 
surfactant into a water molecule and a oil molecule. The 
crucial feature of the Kawakatsu-Kawasaki dynamics, 
which distinguishes it from the Kawasaki spin-exchange 
dynamics for binary alloys, is that the anti-parallel spin 
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model, namely, m = q = 1, p = 0.|— . 

Jan, StaufFer and coUaboratorsEJ extended the model 
further to describe single-chain surfactants with ionic 
heads. According to their formulation, the monomers 
belonging to the ionic heads have Ising spin +2 to mimic 
the presence of electric charge. The repulsive interac- 
tion between a pair of ionic heads is taken into account 
through an (antiferromagnetic) interaction J = — 1 be- 
tween pairs of nearest neighbour sites both of which carry 
spins +2; however, the interaction between all other pairs 
of nearest-neighbour spins is assumed to be J = 1. By 
restricting the range of the repulsive (antiferromagnetic) 
interaction between the "charged" heads to only one 
ttice spacing one is, effectively, assuming very strong 
^^reening of the Coulomb repulsion between ionic heads 
Dy the counterions. 

Starting from an initial state (which depends on the 
ihenomenon under investigation), the system is allowed 
£> evolve following the standard Metropolis algorithm: 
ach of the attempts to move a surfactant takes place 
^rtainly if AE < and with a probability proportional 
exp{-AE/T) if AE > 0, where AE is the change in 
energy that would be caused by the proposed move of 
the surfactant under consideration. 

Gemini with spacer 6 



• Head o Tail/Spacer e 

FIG. 3. Larson-type models of single-chain, double-ch 
gemini surfactants. 
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Jan, StaufFer and coUaboratorsEzI simplified the LVrson 
model by formulating it in terms of Ising-like variabrcs in 
the same spirit in which a large number of simpler Mfetice 
models had been formulated earliertS for the convei^nce 
of calculations. In this reformulation, a classical Ising 
spin variable S is assigned to each lattice site; "J^g^pltation 
(—1) if the z-th lattice site is occupied by a water (oil) 
molecule. If the j-th site is occupied by a monomer be- 
longing to a surfactant then Sj = 1,-1,0 dependi( 
whether the monomer at the jth site belongs to he 
or neutral part. The monomer- monomer interactio, 
taken into account through the interaction betwee 
corresponding pair of Ising spins which is assumed 
non-zero provided the spins are located on the m 
neighbour sites on the lattice. Thus, the Hamiltoni( 
the system is given by the standard form 

H = —J ^ ' SiSj. 

<ij> 

where attractive interaction (analogue of the ferrcBiBBSkling 
netic interaction in Ising magnets) corresponds to J > 
and repulsive interaction (analogue-of antiferromagnetic 
interaction) corresponds to J < OEII. The temperature 
T of the system is measured in the units of J/ks where 
kg is the Boltzmann constant. The Kawakatsu-Kawasaki 
model may be regarded as a special case of the Larson 
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FIG. 4. Different moves of the amphiphiles 



Next, we specify the allowed moves of the surfactants 



6 



for the appropriate sampling of the states of the system 
in a MC simulation. So far as the single-chain surfactants 
are concerned, the moves allowed for the surfactants are 
as follows: (see fig. 4) 

(i) reptation: one of the two ends of each surfactant is 
picked up randomly, with equal probability, and the sur- 
factant is allowed to move forward along its own contour 
by one lattice spacing with the probability mentioned 
above; this move effectively mimics the reptile-like slith- 
ering of the surfactants and hence the name; 

(ii) spontaneous chain buckling: a portion in the middle 
of the tail is randomly picked up and allowed to buckle 
with the probability mentioned above; 

(iii) kink movement: a kink formed by the buckling or 
reptation is allowed to move to a new position with the 
appropriate probability calculated according to the pre- 
scription mentioned above; 

(iv) pull move: this is the reverse of spontaneous chain 
buckling; a buckled part of the tail is pulled so as to make 
it more extended. Each of these moves is possible only if 
the new positions of all the monomers are not occupied 
simultaneously by monomers belonging to other surfac- 
tants. Each surfactant is allowed to try each of the above 
mentioned moves once during each MC step. 

Note that the monomers of the same surfactant as 
well as different surfactants are not allowed to occupy 
the same lattice site simultaneously; this represents a 
hard-core intra-chain as well as inter-chain repulsion 
for monomer separations smaller than one lattice spac- 
ing. Moreover, at any non-vanishing temperature, dur- 
ing the out-of-line thermal fluctuations of the chains, the 
hard-core repulsion leads to steric repulsion between the 
chains. To our knowledge, no potential energies associ- 
ated with the torsion of the surfactant chains have been 
incorporated so far in any work on Larson-type models. 

A microscopic lattice model of double-chain surfac- 
tants (with a single head) in aqueous solution was de- 
veloped by BernardesES by tnqdifying the Larson model 
of single-chain surfactantsEaEZlCj. In terms of the sym- 
bols used above to denote the primary "structure" of 
the microscopic lattice model of single-chain surfactants, 
Bernardes' lattice model of double-chain surfactants, 
with a single hydrophilic head, can be described by 
the symhol-T,nJ\fpHqN'p%n. Very recently, Maiti and 
Chowdhuryc3 have proposed a microscopic lattice model 
of gemini surfactants by extending Bernardes' model so 
as to incorporate two hydrophilic heads connected by a 
spacer . This model of a gemini surfactant can be rep- 
resented by the symbol TmNp'HqSn'HqNpTm where n is 
the number of lattice sites constituting the spacer rep- 
resented by the symbol S and the other symbols have 
the same meaning as in the case of Bernardes' model of 
double-chain surfactants (see fig. 3). For the convenience 
of computation, the Bernardes' model of double-chain 
surfactants as well as Maiti and Chowdhury's model of 
gemini surfactants have also been formulated in terms of 
classical Ising spin variables, generalizing the correspond- 
ing formulation for the single-chain surfactants reported 



in rel0. |_. 

To our knowledge, the work of Bernardes et alxA is the 
only published report of the investigation on the kinetics 
of phase ordering in the Larson model of microemulsions. 
Oil, water and surfactant molecules were distributed ran- 
domly in the initial state of the system. The system 
was allowed to evolve by implementing the moves of the 
surfactants according to the algorithm mentioned ear- 
lier while the molecules of oil and water were allowed to 
exchange their positions according to the Kawasaki-spin 
exchange dynamics mentioned earlier in the context of 
binary mixtures. The time evolution of the mean sizes of 
the clusters of oil, water and amphiphilcs in the system 
were observed up to 200,000 MC steps. 

When the amphiphile concentration (pa = 0.1 and the 
concentrations of oil and water (po = (f>w = 0.45, respec- 
tively, the average length scale < R >ow= [(l/2)(Xo + 
Xw)Y^^ was found to obey the power law < R >ow'^ 
with n ~ 1/3, i.e., the coarsening was found to be gov- 
erned by the same Lifshitz-Slyozov law which governs 
the coarsening in binary mixtures in the absence of am- 
phiphilcs. However, comparing these results with the cor- 
responding results obtained by other groups using vari- 
ous different techniques, we believe that the truly asymp- 
totic logarithmically slow growth regime lies beyond the 
longest time scales of observation in the computer exper- 
iments of Bernardes et al. 

The kinetics of ordering in microemulsions and micel- 
lar solutions containing gemini surfactants have not been 
reported so far. A microscopic model for single-chain 
surfactants at the air-water interface has been devel- 
oped earlier by one-nfjidlj'LJ by appropriately modifjdng 
the Larson modcloEZIc3 of ternary microemulsionsEilEI. 
Later Maiti and Chowdhury replaced the single-chain 
surfactants in the model introduced in ref.a by the model 
gemini surfactants , thereby getting the desired micro- 
scopic model of gemini surfactants at the air-water inter- 
face. A novel entropy-driven phase segregation has been 
observed in computer experiments on a binary mixture 
of chemically identical single-chain surfactants of two dif- 
ferent lengths at the air-water interface. However, its 
dynamical aspects have not been investigated so far. 

C. Hybrid Models 

In this model the binary mixture of A (say, oil) and B 
(say, water) is represented by a continuum field i/'(^) in 
the same fashion as outlined in section 2 in the context 
of the continuum description of binary alloys. However, 
unlike the continuum model discussed in the preceding 
section, the amphiphilcs are treated at the molecular level 
by describing the dynamics of the positions of the center 
of gravity and the orientations of the molecules. That is 
why this model is called "hybrid" . 

In this model the surfactant molecules are modelled 
as "dumbbells" of length I which have two interactions 
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centers at the two ends, one of which is A-philic and the 
other is i?-phihc. Kawakatsu and KawasakicJ assumed 
that the "A-phihc and i?-phihc interactions centers of 
the surfactant have the same chemical species as A and 
B components of the binary mixture, respectively" . Sup- 
pose, the position of the center of gravity of the i-th sur- 
factant molecule is denoted by fi and the unit vector from 
the _B-philic interaction center to the A-philic interaction 
center of the same molecule by the symbol Si. For sim- 
plicity, also assume that Vaa(?0 = ^bb(^ = and 
^b(^) = x{^- Then, the free energy functional is given 
by 



- tps 



(18) 



where 



F^^ = J d'*r[(l/2)Cv,{VVXr)}2 + {-ro)^j^r) + u^\r)] 

(19) 

F^, = Ai.iV, -I- iql/2)J2 J d'^rV-if- rl)s,.VV^(f) (20) 
Fss = J2^2V+{n,) + (/V4)(s-.s-.)VV<^(r,,) 



with 



+ aV4)(44)VVx(r.,)] 



V±(r) =0(r)±x(r) 



where 



P = PA{r) + PBir) 



(21) 

(22) 
(23) 
(24) 

(25) 
(26) 



is assumed to be a constant. The equations of motion 
are given by 



dij{r,t)/dt = L^V'^[dF/6ij] 



dn{t)/dt = ^Lp{dF/dn) 



(27) 



(28) 



dh{t)/dt = -L,[{dF/dsi) - {{dF/dsi).h}h] (29) 

where L's are phenomenological kinetic coefficients. The 
second term on the right hand side of the Eq. (29) arises 
from the constraint \si\ = 1. The Eqs. (27)- (29) were 



solved numerically assuming the forms (t>{r) = — exp(— r) 
and x{'^) = aexp(— r). Equation (27) was solved by the 
cell dynamic method whereas Eqs. (28) and (29) were 
solved by the molecular dynamic method (see the original 
paper for the numerical values of the various parameters). 
Two kinds of initial conditions were used: (a) in case 
of equal volume fractions of the A and B components 
a random bicontinuous structure formed and coarsened 
with time (b) in the case where the volume fraction of the 
B component were three times that of the A component 
a dispersion of droplets was formed which coarsened with 
time. 

In both the situations |-(a) and (b) mentioned above, 
Kawakatsu and Kawasakie3 found a crossover in the tem- 
poral evolution of R{t) from the power law R{t) ~ t^^^ 
to a slower growth. Such a crossover was found to oc- 
cur when the oil-water interface gets saturated by the 
surfactants. j— . 

Kawakatsu et aZ.c3 have also studied phase separa- 
tion in immiscible binary mixtures in the presence of 
surfactant molecules of asymmetric shape using the hy- 
brid model. They observed a spontaneous morphological 
change from a bicontinuous structure to a micellar struc- 
ture during the phase separation process which gives rise 
to the formation of a amphiphilic monolayer at the oil- 
water interface. This is consistent with the known fact 
that the asymmetry of molecular shape gives dsjSj-to spon- 
taneous curvature of amphiphilic monolayerscZrEj. 



D. Purely Phenomenological Models 

This model is an extension of the continuum model of 
binary alloy phase ordering. Naturally, this is an appro- 
priate extension of Eqs. (2)-(4) above. Suppose p{r) is 
the density c^Jjhe surfactant molecules at the location r. 
Laradji et aZ.E^CJ postulated that the effective free energy 
functional for the ternary system under consideration is 
given by 

F[i![r),p{r)] = / d'^r[c{\Itpf - To^^ + wV''' 



^lp ~ + Fs 

(30) 



Here p{r) is the density of the surfactant molecules at 
the location r. c,ro,u^ g, a and /i are phenomenological 
coefficients, is the chemical potential of the amphiphilic 
molecules. is the difference of the chemical potentials 
of water and oil. Note that g is the strength of the cou- 
pling between the two fields p and ip and Eq.(30) is an 
extension of equation (2). The equations of motion in 
this model are given by 



dp/dt = V^{SF/6p)+^pif,t) 



(31) 



(32) 
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where 

< 77(f, t)ri{r',t') = 2kBTV'^5{r- r')5{t - t') (33) 

The above equations (31)-(33) are the generahzations 
of the Eqs. (3) and (4). In the termijBology of the 
Hohenberg-Halperin classification schemeE3 this model 
corresponds to the model D. The surfactant property of 
the amphiphiles is taken into account through 




where s is a phenomenological constant. Laradji et 
aZ.lla solved the Eqs. (30)-(33) numerically (see the orig- 
inal paper for the numerical values of the parameters) 
starting from a random initial condition, i.e., to each grid 
point i^{r) and p{r) were assigned small random values 
around their initial values a,i t — 0. The second mo- 
ment of circularly averaged structure factor, R~^(t), was 
monitored as a function of time t. The main results of 
their investigations are as follows: (i) The location of the 
peak of the structure factor moves initially to small k as 
time passes, thereby indicating coarsening. However, the 
coarsening seems to come to a halt at very late stages be- 
cause the peak position was observed to become static at 
a fixed k = fee ^ 0. Moreover, the larger is the concen- 
tration of the surfactants, the smaller is the final size of 
the oil- rich (or water- rich) domains, (ii) R{t) oc {logt)^, 
although estimation of y was not carried out by these 
authors. The slow growth observed in this study is con- 
sistent with the corresponding results of the microscopic 
models in the preceding section—-. 

Recently, Ahluwalia and PuriE3 have made an attempt 
to study the kinetics of phase ordering in ternary mi- 
croemulsions using a CDS approach. The results are 
consistent with those obtained by the other methods dis- 
cussed so far. 



IV. SUMMARY AND CONCLUSION 

For the description of those phenomena which are 
strongly dependent on steric interactions, the Larson 
model is believed to be more realistic than the lattice 
model of Kawakatsu and Kawasaki. However, if these 
two models belong to the same dynamic universality 
class, then the long-time regime, where Kawakatsu and 
Kawasaki have observed non-algebraic slow growth, must 
be-iying beyond the longest runs made by Bernardes et 
al.a 

It is worth mentioning here that the hybrid model 
takes into account the surfactant property of amphiphilic 
molecules in a much more realistic manner than that in 
the purely phenomenological models. Nevertheless, the 
hybrid model suffers from the shortcomings that the ex- 
cluded volume of the surfactants is not taken into account 
because of the assumption that pA{r)+pB{'r) = constant. 



The amphiphilic molecules in ternary microemulsions 
and micellar solutions may be regarded as impurities 
added to an immiscible binary mixture of oil and water. 
But, these are not "quenched" impurities and, therefore, 
the process of phase ordering is not driven by thermally 
activated motion of the oil-water interface. Almost all 
the numerical works published so far demonstrate that 
during the late stages of the coarsening process the oil- 
water interface is saturated by the surfactants leading to 
a reduction of the oil- water interfacial tension which gives 
rise to the crossover from the Lifshitz-Slyozov-like power- 
law growth to a slower growth of the ordered regions. 
This is consistent with the results obtained by incorpo- 
rating the effects of the surfactants in the Lifshitz-Slyozov 
theory for the [kinetics of phase ordering in immiscible bi- 
nary mixturesE3. 

Experimental study of the kinetics of ordering in 
ternary microemulsions and micellar solutions is difficult 
because this phenomenon is too fast to observe over time 
scales which would be sufficient to extract the growth 
laws. A — B diblock copolymer in a A/B binary ho- 
mopolymer blend behaves effectively like a surfactant be- 
cause the A [B) subchain of the copolymer dissolves pref- 
erentially in the A-rich (B-rich) phase Ell. Since the phase 
separation process in polymer systems is much slower 
because of entanglements of the polymer chains, the ki- 
netics of ordering in a binary polymer blend containing 
an amphiphilic block copolymer can be utilized for ex- 
perimental investigations. The results of such investiga- 
tions can be compared with the theories of the kinetics 
of phase ordering in ternary microemulsions and micellar 
solutions. 
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